epi draft 



Interaction of two-dimensional plasma crystals with upstream 
charged particles 

C.-R. Du, V. NosENKO, S. Zhdanov, H. M. Thomas, and G. E. Morfill 

Max Planck Institute for Extraterrestrial Physics, Giessenhachstrafie, 85748 Garching, Germany 



o 

Oh 



PACS 52.27.Lw - Dusty or complex plasmas; plasma crystals 
PACS 52 . 25 . Fi - Transport properties 
PACS 52 . 27 . Gr - Strongly-coupled plasmas 

Abstract. - Two-dimensional plasma crystals are characterized by a strong up-and-down asym- 
metry not only due to gravity but also due to the presence of plasma flow at the location of 
particles. We study for the first time the interaction of a single-layer plasma crystal with charged 
extra particles located above it (upstream of the flow of ions). Upstream extra particles tend 
to move between the rows of particles in the crystal, accelerate to supersonic speeds, and excite 
attraction-dominated Mach cones and wakes in the crystal. 



Qh Introduction. — A complex plasma is a weakly ion- 
c/^ized gas containing ions, electrons, as well as small solid 
.^particles The particles are usually negatively 

^charged due to the high thermal speed of electrons. The 
^^harge of particles ranges from a few hundreds to several 
O thousands of elementary charges depending on the parti- 
1— ^le size and discharge conditions. In ground-based exper- 
^ iments the particles can levitate in the pre-electrode area 
^against gravity by a strong electric field. Under certain 
experimental conditions they can be confined in a single 
Cn| layer and self-organize in a triangular lattice with hexago- 
nal symmetry. Such system is known as 2D plasma crystal. 
' ^The in-plane interaction between particles in this system 
r is well described by screened Coulomb (Yukawa) poten- 
^^tial with the screening length defined mainly by electrons 

• • Since the discovery of plasma crystals [6||8], various 
, ^ experiments have been performed with 2D plasma crys- 
^^tals including melting [9|p!Q], recrystallization 11 , defect 

^transport 12 , etc. In the course of such experiments. 



^ one often observes extra particles apart from the main 2D 
lattice layer after injecting particles into plasma. Those 
particles, which can be agglomerates or contaminations, 
sometimes move at a high speed, disturb the lattice, and 
can easily spoil the desired experiments by creating wave 
patterns within the lattice layer. When the extra parti- 
cle speed is higher than the sound speed of the lattice, 
the disturbance forms a Mach cone. The Mach cones and 
wakes associated with extra particles moving beneath the 
lattice layer were well studied in the past decade 13-17 . 



This phenomenon can be used for diagnostic purpose. In 
fact, by measuring the angle of a Mach cone and the speed 
of the source of disturbance, the sound speed of the crys- 
talline lattice can be directly estimated. This method was 



used for 2D as well as 3D plasma crystals 17 -19 . Besides 
measuring the sound speed of a plasma crystal, one can 
also use the extra particles to heat the crystalline lattice 

However, in many cases one needs an undisturbed 2D 
plasma crystal to perform some delicate experiments, e.g., 
to investigate the dynamics of a perfect crystal. It is rel- 
atively easy to get rid of the extra particles beneath the 
lattice layer (downstream of the flow of ions). In practice 
one can drop those extra particles on the bottom electrode 
by reducing the discharge power at higher pressure. We 
call this process "purification" . 

In the experiments performed in our laboratory, Mach 
cones and related wakes were sometimes observed in 2D 
plasma crystals even after purification. Such wave pat- 
terns happened to be induced by extra particles located 
above the lattice layer (upstream of the flow of ions), and 
showed many different features. 

In this Letter we report for the first time on the obser- 
vation of channeling and leapfrog motion of upstream ex- 
tra particles, accompanied by the excitation of attraction- 
dominated wakes in the lattice. 

Experimental setup and conditions. — The exper- 
iments were performed in a modified Gaseous Electronics 
Conference (GEC) rf reference cell 21 , see Fig. [l] Ar- 
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Fig. 1: Sketch of experimental setup with a modified GEC 
chamber. The bottom electrode is powered by rf generator at 
13.56 MHz, the upper grounded ring and the chamber walls 
(not shown here) serve as the counter-electrode. Laser 1 shines 
a horizontal laser sheet to illuminate the lattice layer of 2D 
plasma crystal while laser 2 illuminates an extra particle out- 
side the lattice layer. The optical axes of the lasers are actually 
oriented at 90° to each other. Particle motion in the lattice 
layer as well as of the extra particle can be simultaneously 
recorded by a video camera with top view. 



gon plasma was sustained using a capacitively coupled rf 
discharge at 13.56 MHz. The input power was set at 20 W. 

We used monodisperse melamine formaldehyde (MF) 
and polystyrene (PS) particles to create 2D plasma crys- 
tals suspended above the bottom rf electrode. The MF 
particles have a diameter of 9.19 ±0.09 /xm and mass den- 
sity of 1.51 g/cm^, while the PS particles have a diame- 
ter of 11.36 ± 0.12 /im and mass density of 1.05 g/cm^. 
Gas pressure was maintained at about 0.65 Pa; the corre- 
sponding neutral gas damping rate was 7 0.77 s~^ for 
MF particles and 7 0.91 s~^ for PS particles 22 . Fur- 
ther experimental parameters are listed in Table [l] The 
lattice layer was illuminated by a horizontal laser sheet 
shining through a side window of the chamber. A high- 
resolution video camera (Photron FASTCAM 1024 PGI) 
was mounted above the chamber, capturing a top view 
with a size of 42 x 42 mm^, as sketched in Fig. [l] The 
recording rate was set at 60 frames per second. 

Mach cones and wakes. — We performed two sep- 
arate experiments intended to observe two different types 
of Mach cones. For the sake of simplicity, in this Letter 
we name the Mach cones in the lattice observed after its 
purification the "type I" Mach cones. In contrast, the well- 
known Mach cones excited by extra particles beneath the 
lattice layer are called the "type 11" . To observe the type II 
Mach cones, we did not purify the plasma crystal so that 
big extra particles remained in the discharge. To make 
the results comparable, two cases with the similar Mach 
number^ but related to two different excitation sources 
were selected, see Fig. |2j The particle areal number den- 



^ The Mach number is defined as M = v^jCi^^ where is 
the speed of disturbance (extra particle) and Cj^ is the longitudi- 
nal sound speed of the crystalline lattice. 



15 



10 



' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 
(a) 



n 

-(mm ) 
4.0 5.0 



1 1 1 1 I 1 1 1 I I 1 1 1 I I 1 1 1 I I 1 1 1 I I I 



15 



' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I 



(b) 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 



1 2 3 4 5 n(mm" ) 



8 wT'irr- 



7 

6 

>^ 4 
3 
2 



-T^i^v^^y^^M"^ (C) 

^ ^ f ff — ^ - 

V V V ^ fX-^ < K r 

A - ^ r 



^ y 
^ y 



t 



I i>i I n^i I i>i I i-H I I I lyri ri^l I I 



7 8 9 
X (mm) 



10 




net force 50 fN — > 



50 fN^ 



Fig. 2: Comparison of the particle areal number density n and 
net force between (a), (c) type I (experiment 1 in Table [T]) and 
(b), (d) type II (experiment 2 in Table [T]) Mach cones. The 
maps for type I and type II Mach cones are averaged from 
data for 30 and 20 consecutive video frames, respectively. The 
insets show the density profile (solid lines) along the symmetry 
axis (marked by the dash-dotted lines). The dashed lines in the 
insets represent the density profile for the undisturbed lattice. 
The net force Fn = ma -h 7mv acting on the particles in the 
lattice at the apex of the Mach cone [marked by the dashed 
rectangle in (a) and (b)] is shown in panels (c) and (d). The 
orange circles highlight the apex of the Mach cones. 



sity map and the net force vector field plot corresponding 
to the type I Mach cone are shown in the left panels in 
Fig. |2j The wings of the Mach cone are clearly recogniz- 
able in both plots. The half-opening angle between the 
wings ("cone angle" in what follows) is measured to be 
fi = 55° ± 10°. By using the well-known Mach cone rela- 
tion sin /i — M~^, where M > 1, one can obtain M 1.2. 
The measured Mach number agrees well with the value 

Note 



estimated using the phonon spectrum method 23 



that the lateral wakes are also well resolved in the density 
map; these can be used for diagnostic purpose [24]. 

In the density map. Fig. |2|a), there is a density in- 
crease at the apex of the type I Mach cone. Comparing 
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Table 1: Experimental parameters" including the particle material and diameter d, charge number Z — —Q/e, interparticle 
distance^ A, screening length Ad, longitudinal sound speed Cl, extra particle speed Vd, Mach number M, and Mach cone type 
(see text for details). The particles were suspended in an argon discharge at the pressure of 0.65 Pa and discharge power of 



20 W, 



experiment 


particle 


d 


Z 


A 


Xd 


Cl 


extra 




M 


Mach cone 




material 


(/im) 




(/im) 


(/im) 


(mm/s) 


particle 
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type 
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MF 


9.19 


15400 


520 


390 


27 


upstream 


32 


1.2 


I 
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MF 


9.19 


16300 


560 


600 
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downstream 


42 


1.3 


II 
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MF 


9.19 


15000 


610 


300 


17 


upstream 


24 


1.4 
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4 


MF 


9.19 


16100 


560 


400 


26 


downstream 


42 


1.6 


II 


5 


PS 


11.36 


19000 


650 


560 


27 


upstream 


29 


1.1 


I 



" Error bars: ±13% for Z, ±45% for A^, ±15% for Cl, ±5% for Vd, and ±15% for M. ^ A is obtained from the first peak of 
the pair correlation function. The crystalline lattice is slightly inhomogeneous (2% at the most). 



to the undisturbed valuq^ the particle density increases 
from 4 mm~^ to 4.5 mm~^. The reason for that is simple: 
The particles in the lattice layer are dragged toward each 
other, resulting in an increase of the local particle density. 
This is clearly seen in the net force vector field shown in 
Fig-ic). 

As to the "conventional" type II Mach cone, the density 
at the apex drops down sharply from 3.7 mm~^ (undis- 
turbed particle density) to 1.9 mm~^, producing a hole in 
the lattice clearly seen in Fig.[2jb). The density drops be- 
cause the lattice particles at the apex of the Mach cone are 
repelled by the extra particle, see Fig.|2|d). The repulsion 
is caused by the Yukawa repulsive force between the extra 
particle and the particles in the lattice layer, all of which 
are negatively charged. The cone anglais /i = 50° ± 6°. 

Extra particle diagnostics. — In order to trace the 
extra particle originating the Mach cone, we installed a 
second laser at 90° to the first one. This laser shines a laser 
sheet, which is parallel to the first laser sheet, as shown 
in Fig. [l] The height of the laser can be finely adjusted 
with accuracy of 10 /im. Using two lasers simultaneously 
and registering the light scattered by the particles with the 
same camera allowed us to image the lattice layer together 
with the extra particle. (The depth of field of the camera 
lens is larger than the distance between the extra particle 
and the lattice layer.) The trajectory of the extra particle 
is therefore overlapped with the trajectories of particles in 
the lattice layer, see Figs. |3|a), (b). That fits well with 
our goals. 

The extra particles beneath the lattice layer moved at 
much lower heights. They could be observed if we lowered 
the height of the second laser. For brevity, we call these 
particles "downstream" because they are located down- 
stream of the ion flow in the (pre) sheath comparing to the 
particles in the lattice layer. As shown in Fig. [sj^b), the 
trajectory of this extra particle beneath the lattice layer is 
well resolved as well as the displacement of individual par- 

^ Note that the unperturbed crystalline lattice inhomogeneity is 
small compared to the density variation caused by the extra particle. 
^Measured in the particle speed map (not shown here). 



tides in the lattice layer. The extra particle trajectory is 
rather smooth and apparently not influenced by the local 
structure of the lattice layer above it. By measuring the 
height difference between two laser sheets, we estimate the 
vertical distance of extra particle from the lattice layer to 
be Ah 800 jam. The Mach cone can already be seen in 
Fig. [sj^b). However, the cone structure is better resolved 
in the velocity vector field. Fig. |3jd). It is clear that the 
extra particle repels the particles in the lattice layer right 
above it due to the Yukawa repulsion. The apex of the 
Mach cone follows this extra particle closely. 

It is more difficult to find the source of type I Mach cone. 
We scan the second laser vertically with small height steps 
to catch the trace of any possible source of disturbance. 
It turns out that the source of type I Mach cone is an ex- 
tra particle moving above the lattice layer. For the reason 
explained above, we call these particles "upstream". The 
relative height of these particles above the lattice layer is 
estimated to be Ah 200 /im. There are several reasons 
why it is hard to visualize or in fact even notice the ex- 
istence of such particles: (i) This type of extra particles 
appears in a plasma crystal much less often than down- 
stream particles, (ii) these particles are generally smaller 
in size, resulting in a very dim trajectory recorded by the 
video camera, (iii) the influence of such particle on the 
lattice layer is relatively small. As we already see in the 
density maps in Fig. [2j the density variation caused by 
an upstream particle is much smaller then that caused by 
downstream particle. Thus it is easy to overlook such ef- 
fect without careful analysis. 

As a rule, the upstream particles appear to move slower 
than downstream particles. Their average velocity de- 
pends on the gas pressure, discharge power, particle size, 
etc. Such dependence is not the focus of this Letter and 
will be reported elsewhere. 

Extra particle — lattice interaction at the cone 
apex. — The computed velocity of particles in the lat- 
tice layer with Mach cones excited by upstream and down- 
stream particles are shown in Fig. ^c) (type I) and 
Fig. [sj^d) (type II), respectively. It is immediately clear 
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Fig. 3: Trajectories of upstream (a) and downstream (b) extra 
particles (experiments 3 and 4 in Table [T]) . The particle po- 
sitions are presented by superimposing a series of consecutive 
experimental images for 0.45 s (a) and 0.22 s (b) (color-coded 
from blue to red). The trajectory of the upstream particle, 
showing a clear zig-zag feature, is magnified in the inset in (a). 
The scattering angle is a = 25° di 3°. Panels (c) and (d) show 
snapshots of positions (black dots) and velocities (black scaled 
arrows) of particles in the lattice layer in a single frame as well 
as positions (red dots) of upstream and downstream particles 
and their velocities (red arrows, scaled 1:5), respectively. In 
(c) and (d) the Mach cone is highlighted by orange lines. The 
cone angles are 44° ± 12° (c) and 38° ± 8° (d). 



that they are strikingly different. First, in the vicinity of 
the apex of type I cone, the lattice layer particles move 
toward the upstream particle [Fig. |3jc)], whereas they are 
visibly repelled by the downstream particle [Fig. [sj^d)]. 
Second, the apex of the Mach cone in Fig. [sjc) has a elon- 
gated shape along the direction of the motion of the extra 
particle. The conjunction of the wings of the cone is there- 
fore located slightly behind the moving upstream particle 
at a distance ~ (2 — 3)A. In contrast, the apex of the 
type II cone follows the extra particle position closely, as 
shown in Fig. |3|d). 

Importantly, the force field in Figs. [2|c),(d) and the 
velocity field in Figs. |3jc),(d) are similarly "polarized", 
i.e., both fields are simultaneously directed toward or 
away from the instantaneous location of the extra par- 
ticle. Based on what we observed, the extra particle - 
lattice layer interaction is dominantly attractive for type I 
cone (upstream particle) or dominantly repulsive for type 
II cone (downstream particle). 

We attribute the particle-lattice layer attraction to the 



ion wake formed underneath the upstream extra particl^ 
Positive ions that concentrate locally below the extra par- 
ticle exert an attractive force on the negatively charged 
particles of the lattice layer. This phenomenon is well 
known as the ion-wake effect causing, e.g. the particles 



to pair |2,25-28j. The particle-lattice repulsion resulting 
from the Yukawa repulsive force at the apex of the type 
II cone is a well-established effect 
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Finally, the situation appears to be more complicated, 
and the particle - lattice interaction cannot be reduced 
only to the dominantly-attractive or dominantly-repulsive 
interactions. The examples discussed below can be prop- 
erly addressed only suggesting the presence of competing 
repulsive-attractive interactions (which are also considered 
to be important for colloids [29].) 

Channeling effect. — There is one more surprising 
observation regarding upstream particles: The long-term 
behavior of these particles demonstrates anomalous trans- 
port properties and elements of "strange kinetics" 
For instance, the particle is apparently able to channel 
between two aligned rows of particles in the crystalline 
layer ("wall particles" below), as seen in Fig. |3ja). 

Confinement force. The upstream particle is well con- 
fined in the channel. Comparing to the smooth trajectory 
of the downstream particle. Fig. [sj^b), the trajectory in 
Fig. [3]; a) has a zig-zag shape [see inset where a part of the 
trajectory in question is magnified]. In other words, the 
extra particle, interacting with the crystal particles com- 
prising the channel, evidently bounces and hence the effect 
of the wall particles on the extra particle is dominated by 
repulsion. By measuring the scattering angle a, one can 
estimate the confinement force Fc by using the relation 
FcAt = m < Vd > {I — cos a), implying momentum con- 
servation, where At o:^ 0.08 s is the time of scattering, m 
and < Vd > are the mass and mean (longitudinal) veloc- 
ity of the extra particle, respectively. In the experiment 
presented here, the angle is measured to be a = 25° ± 3°, 
Vd = 29 ^2 mm/s, and the confinement force is estimated 
as Fc ^ 21 fN. This value is about twice the neutral gas 
friction force — 13 fN. 

Non-reciprocal interaction. The extra particle, pass- 
ing through the channel, in turn exerts the force on the 
wall particles. This force deforms the lattice cells when the 
extra particle moves through. Studying this deformation, 
a certain conclusion can be made about the wall particle - 
extra particle interaction. From Fig. (sja) it follows imme- 
diately that this interaction is non-reciprocal because the 
wall particles behave as if they are attracted to rather than 
repelled by the similarly charged extra-particle. This kind 
of non-reciprocity is indeed easy to address taking into 
account the ion wake (ion focus [2^) formed beneath the 
upstream particle. 



^ As is the case for any particle levitating in a discharge 
(pre)sheath. 

^ Channeling is a process that constrains the path of a charged 
particle in a crystalline solid (31]. 
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Fig. 4: Channeling of an upstream particle (a) and related 
cell deformation (b) (experiment 5 in Table [TJ. Panel (a): 
black dots represent the lattice particles, colored dots represent 
positions of the extra particle and the "wall" particles color- 
coded from blue to red for 0.66 s. Panel (b), inset: A magnified 
single cell [marked by a dashed-line rectangle in (a)] with color- 
coded particle positions (from blue to red with a time step of 
4 ms). The dashed line is the least squares linear fit. 

Lattice cell deformation. In order to investigate the 
deformation of the lattice cell caused by the passage of 
extra particle in more detail, we performed another exper- 
iment with higher recording rate of 250 frames per second 
(experiment 5 in Table [T]). Accordingly, we used particles 
with larger diameter so that they can be well illuminated 
with the same illumination laser power without a signif- 
icant change of the recording quality. An (unexpected) 
advantage of this case is that the extra particle above the 
lattice layer can even be traced using the same laser. This 
implies that the height difference is roughly of the laser 
sheet width or even smaller. Ah < 100 /im. This height 
difference is about one-fifth of the characteristic length 
Le = E/E'h ^he levitation height in the (pre) sheath, 
where E is the local electric field The expected dif- 
ference in particle mass Am/m ~ 1.5Ah/LE is less than 
30-35%. 

The extra particle track and the accompanying cell de- 

^ Note that both PS and MP particles levitated at approximately 
the same height. 



Fig. 5: Interaction of upstream particles with defects in the 
crystalline lattice. The particle positions obtained for the time 
interval of 0.4 s are shown as the solid dots color-coded from 
blue to red. Initially in both cases the extra particles move 
in the channel (highlighted by two parallel solid lines). The 
initial velocities are approximately the same = 26di 1 mm/s. 

(a) - single point defect scattering; the scattered particle is 
captured in a new channel (highlighted by two parallel dashed 
lines) where it moves with the velocity Vd — 27 ± 1 mm/s. 

(b) - interaction with a large-angle grain boundary; the extra 
particle, being initially scattered by a single defect, left the 
channel and penetrated the chain of defects, jumping from one 
position to another with the mean velocity Vd — 19 ^3 mm/s. 
The 7-fold cells are marked by triangles, and 5-fold cells are 
marked by squares. 

formation are shown in Fig. [4j The extra particle moves 
along the channel, which even bends slightly at x > 14 mm 
as shown in Fig.[4]^a). The deviation from the major track 
is small and a zig-zag feature is barely visible in this figure, 
indicating stronger confinement as compared to Fig. [sja) . 
Let us focus on the dynamics of a single cell, see the in- 
set in Fig. [4]^b). When the extra particle enters (through 
the right boundary in the case considered) the cell de- 
forms as indicated by the color-coding in the inset. In 
order to quantify this deformation, we introduce the rel- 
ative cell deformation (A2 — A^)/ where A2(^d) is the 
undisturbed (disturbed) length of the cell side transversal 
to the extra particle track; the relative position of the ex- 
tra particle inside the cell S/ Ai^ where Ai is the length of 
the side of the undisturbed cell along the track, and S is 
the displacement of the extra particle measured from the 



right boundary of the ceh, as shown in Fig. [4j The defor- 
mation depends hnearly on the relative position of extra 
particle: 8//\i = k{A2 - Arf)/A2, where k = 0.14 ± 0.01. 
Note also that in fact the deformation does not only de- 
pend on the relative position of the extra particle, but also 
on its velocity: At a fixed (5/Ai, the slowing-down of the 
particle leads to a higher deformation of the lattice up to 
30%, as seen in Fig.|4];b). 

Interaction with lattice defects. The character of the 
upstream extra particle motion depends on the local struc- 
ture of the lattice layer. Sometimes the extra particles 
move in irregular trajectories as we observed in a number 
of experiments (two examples are shown in Fig. |5|. This 
can be explained by the channel distortion. For instance, 
when the channel is blocked by a structural defect the 
extra particle collides with the defect and scatters leav- 
ing the channel. Depending on the relative kinetic energy 
the collision might be elastic or inelastic. After the colli- 
sion the particle is often again involved in the channeling 
process, see Fig. [sja), being again captured by and accel- 
erated along a channel with a new orientation. Sometimes 
it is a challenge for an extra particle to find a new suit- 
able channel, and it remains quasi-free "leapfrog" jumping 
for a longer time from one position to another one but at 
an essentially lower velocity, see Fig.jsjb). This requires 
further careful analysis. 

To conclude, the upstream extra particles happened to 
be extremely useful tool to effectively test the anomalous 
kinetics effects and the particle-lattice interaction. Those 
particles appeared to be weakly, quasi-elastically interact- 
ing with a strongly coupled 2D complex plasma. That 
allowed us to explore the anomalous long-term channel- 
ing, leapfrog motion, attraction-dominated wakes in the 
lattice, all at the "atomistic" level. 

* * * 
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